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The kinematics of the decay of a bound proton is governed by the proton spectral function. 
We evaluate this quantity in ^^O using the information from nuclear physics experiments. It also 
includes a correlated part. The reliability of this evaluation is sufficient to open the possibility of 
correlated cuts in the missing mass and momentum variables in order to identify the decay events 
from the bound protons with a possible increase of the signal to noise ratio. 



I. INTRODUCTION 



Proton decay is an important window for theories be- 
yond the standard model. Several decay channels have 
been experimentally explored, leading to lower limits for 
the corresponding proton lifetimes. One of them is the 
pionic decay mode into a neutral pion and a positron. 
The signature of a decay process is that the sum of the 
four momenta of the decay products should reconstitute 
the proton four-momentum. For a free proton at rest this 
means a vanishing total three-momentum and a missing 
mass equal to the free proton mass. The pionic decay 
experiment at Kamiokande [ll, Q is performed with wa- 
ter Cerenkov detectors in such a way that 8 out of 10 
of the decaying protons are embedded in an oxygen nu- 
cleus, where their spectrum is modified. In a pure shell 
model description corrections are applied to incorporate 
the shell model momentum distribution and correspond- 
ing binding of the shell model orbits of the oxygen nu- 
cleus. However, Yamazaki and Akaishi Q pointed out 
that this procedure does not take into account the corre- 
lations of the decaying proton with the neighboring nu- 
cleons. Using a correlation function deduced from the 
Reid soft core potential, they evaluated the effect on the 
invariant mass spectrum which acquires a broad low en- 
ergy tail representing ~ 10% of the total decay. 

It is customary in the decay problem to introduce 
the missing three momentum Pmiss = '^iPi, sum of the 
momenta of the decay particles, and the missing mass, 
^miss = {^iEiY — {^iPiY- These quantities also refer 
to the decaying proton. They are related to the mo- 
mentum and energy of the residual nucleus which can 
be in an excited state defined by: Pmiss = Pa-i and 
YiiEi — E\_^. Since the state of the residual nucleus 
for each decay event is not known, there is no strict con- 
straint to identify a proton decay. In the pure shell model 



case where the smearing is already present, the effect is 
rather mild and controllable. But as shown in Ref. Q 
this is not the complete story. It is the aim of this work to 
evaluate these distributions. The issue at stake is if the 
broadening is so big that a substantial fraction of genuine 
decay event are lost in the background. In other words 
the lower limit on the proton lifetime deduced from the 
absence of events within a certain domain in the missing 
momentum and missing energy variables must take into 
account the portion of decay events outside this domain. 
It is therefore useful to have the best possible probability 
distributions. Our result applies to any decay channel 
but for illustrative purposes we will often refer to the pi- 
onic channel. It is valid as well for neutron decay with 
disappearance of the hadron. 

In order to give a feeling for the importance of the mod- 
ification introduced by correlations, in terms of particle- 
hole (ph) excitations, correlations translate into the ex- 
istence of 2p2h excited states mixed into the nuclear 
ground state. The decay of a correlated particle leaves 
the nucleus into an excited state with one hole and one 
particle-one hole, for which we want to evaluate the ex- 
citation energy. Beyond the energy associated with one 
hole creation as evaluated in the shell model, there is the 
energy of the particle- hole which can be approximated by 
Eph = Pph/'^^ where Pph is the momentum exchanged 
between the correlated pair. Neglecting the momentum 
of the hole which has a relatively narrow distribution, 
Pph is also the opposite of the missing momentum. 

The missing mass square is then given by 
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where M* is the nucleon mass reduced by the energy 
necessary for the hole creation. The missing mass square 



evolves approximately parabolically with P„ 
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lustrate the expected effect let us take an approximate 
value M* ~ 900 MeV. For a typical exchanged momen- 
tum Pph = 300 MeV/c, which is also the value of the 
missing momentum, the missing mass value turns out to 
be 800 MeV. These two missing values happen to be on 
the border line of the domain in which Refs. [U, [3] inter- 
pret an event as a proton decay one (no event in fact fell 
into this domain) . We therefore expect that in the analy- 
sis a substantial fraction of the correlated decays escapes 
detection. Moreover, the future experiments aiming to 
improve the current limits on the proton decay will have 
to introduce even tighter cuts to avoid the background 
due to the atmospheric neutrino interactions making the 
effects discussed even more important. The above quali- 
tative argument is made quantitative in the next section. 

II. PROTON SPECTRAL FUNCTION 

Since the decay of the bound proton occurs instanta- 
neously on the scale of nuclear interactions one can ex- 
press the quantities relevant for the bound proton decay 
in terms of the nuclear spectral function 

SA{k,E) = \{i;A\a{k)5{HA-i - E)a\k)\^A)\^ (2) 

which describes the probability to find a nucleon in the 
nucleus with momentum k and produce a residual A — 1 
system with excitation energy E after an instantaneous 
removal of this nucleon. The spectral function is related 
to the single nucleon momentum distribution as 

J dESA{k,E) ^ UAik), (3) 

and it is normalized as 

J dEdkSA{k,E) = 1. (4) 

In order to resolve the spectral function at the high res- 
olution relevant for the proton decay one needs to use 
probes which transfer large energies and momenta, above 
1 GeV, to the nucleons in the nuclei. Such studies were 
performed in the last few years using proton and electron 
beams of high energies. 

It was observed Q that the ratios of (e, e') cross sec- 
tions off nuclei and the deuteron i^He) are independent 
of X, for 1.3 < a; < 2 and > 1.5 GeV^ correspond- 
ing to the kinematics where the electron can scatter only 
off the correlated nucleon - nucleon pair with internal 
momenta > 300 MeV/c. Moreover in {e,e'p) or {p,2p) 
reactions on nuclei at large Q^, a strong correlation was 
observed between the emission of a fast proton and that 
of a nucleon (predominantly neutron) in the opposite di- 
rection Q . These studies confirmed theoretical expec- 
tations of the presence of significant short-range correla- 
tions (SRC) in nuclei - for instance in ■^'^C the probability 
Pi2Q to find a nucleon with momentum > 300 MeV/c is a 



factor of ~ 5 ± 0.5 larger than in the deuteron. The cur- 
rent models of the deuteron give Pd in the range 3 -;- 4%, 
and this corresponds to Pi6q = .15 .2. The data also 
support the expectation that most of this probability is 
due to the pn - tensor correlations (see e.g. 0,[1|), which 
are a specific case of 2p-2h excitations. For a review and 
detailed references see [§]. 

In the many-body models of nuclei with realistic NN 
potential the high momentum component with momenta 
between 300 and 600 MeV/c originates from the interplay 
of attraction and repulsion at distances < 1.2 fm. Hence 
we have used two spectral functions [l0|, Ull calculated in 
such models to analyze the effect of the nuclear structure 
on the detection of the bound proton decay. 

For the purposes of the analysis of the proton decay 
events it is convenient to choose as variables the three- 
momentum of the decaying proton and the square of 
the bound proton mass V — M^^^^ which fixes its off- 
shellness. The two spectral functions of Ref. [lol,[Tl| have 
an uncorrelated part, So{Pm,iss, E), and a correlated one, 
Si{Pmiss, E), where E is the proton removal energy de- 
fined as £' = Ea — E'^_^. In the first model we use Sq 
as calculated with the Skyrme force and renormalized by 
a factor 0.8. The correlated part, Si represents 20% of 
the total spectral function and it is given by the model 
of 10]. In this model / d^kdESi{k, E)0{k ~ fco) = 11% 
where fcp = 300 MeV/c; (the second model we consid- 
ered provides the same result as far as this quantity is 
concerned). The ratio of ^^O and deuteron high momen- 
tum components in these models varies in the range of 
3-^6 for 300 < k < QOOMeV/ c which is rather close to 
the value of the ratio ~ 5 obtained from the analysis of 
the hard phenomena and in particular x > 1 data (see 
the review in Q ) . Smaller value of the total probability 
than in a phenomenological estimate is mainly due to a 
later onset of the dominance of the short-range correla- 
tion regime. The model of Ref. [l^ for the correlated 
part of the spectral function is based on the notion of 
the factorization of the two-body momentum distribu- 
tion for high values of the relative and small values of 
the center of mass momenta of the pair and it is valid 
in this regions; this factorization was justified within a 
many-body approach in Ref. [T^l and shown to hold for 
^^O within the many-body calculation of Ref. It 
also gives a correct dependence for the center of mass of 
the correlated pair, as measured in Refs. ^5^ and [6^. We 
choose for the relative motion of the pair in the two-body 
momentum distribution a parametrization which repro- 
duces well the high momentum tail of the deuteron in the 
region of interest and leads to a good description of the 
high momentum tail of ni6o{k) (l4|. 

The spectral function of Ref. [il| has been obtained 
within the Local Density Approximation [l5| . in which 
the (e, e'p) data on single nucleon knock-out at low miss- 
ing energy [16] is combined with the results of accurate 
theoretical calculations of the nuclear matter spectral 
function at different densities [l3]. A direct measure- 
ment of the correlation component of the spectral func- 
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FIG. 1: Top: the uncorrelated proton spectral function of 
^^O, Pmiss So{Pmiss,V), as Calculated with the Skyrme force 
and plotted vs. missing momentum Pmiss and virtuality, V — 



Bottom: the correlated proton spectral function of 



= from llOj, Pii,,Sl{Prmss,V). 



tion of ^^C, obtained measuring the (e, e'p) cross section 
at missing momentum and energy up to ^ 800 MeV and 
~ 200 MeV, respectively, has been recently carried out 
at Jefferson Lab by the E97-006 Collaboration The 
data resulting from the preliminary analysis appear to be 
consistent with the theoretical predictions based on the 
spectral function of Ref. [ll| . 

The quantity Fj^i^s'^'o (-Pmiss, V) is shown in Figs. [HH 



in three dimensional plots, and P^j^ggSo{Pi 



is 



shown in Figs. [31l3]in contour plots, for the two consid- 
ered models. The strength is concentrated over three or 
four stripes in the Pmiss, M^iss plane in the two cases; 
they correspond to the occupied shells of ^^O: the P1/2, 
P3/2 and 5*1/2 states in the case of calculation with the 
Skyrme force and to an additional P3/2 state in the case 
of Ref. [m. 

The energies and widths of the occupied states in the 
two models we considered are the following. In the 
first model, as calculated from the Hartree-Fock Skyrme 
model with shell model parameters which describe the 
(p, 2p) and {p,pn) data of Refs. they are 12.06 

MeV with a width ~ 5 MeV for the P1/2 state, the P3/2 
state has 18.63 MeV with a width of 5 MeV; the S state 
is quite broad with a width of ~ 40 MeV for an energy 



FIG. 2: As in Fig. [T] using the model of Ref. 
the uncorelated and correlated parts. 



for both 



37.96 MeV. In the second model, Ref. [uj, we have the 
following energies, with roughly the same widths: P1/2 
at 12.50 MeV, P3/2 at 18.75 MeV, a second P3/2 at 23.00 
MeV and 51/2 at 42.50 MeV. The model takes into ac- 
count a rather small contribution of higher energy exci- 
tations where intermediate states with one particle in the 
continuum and one hole in the final state. 

The behavior in the momentum Pmiss can be inferred 
from the expression of V which in the uncorrelated case 
is: 
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Expanding the square in Eq. ([5]) and neglecting the P, 
term we obtain: 
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where are the values of the proton shells energies. The 
coefficients Ca are 1.07 and 1.06 for the P and S pro- 
ton shells, respectively. As for the correlated spectral 
function, P^^^^Si{pmiss,V) is represented in Figs. [H [5] 
in three dimensional plots, and P^j^^ggSi{pmiss, Mmiss) is 
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FIG. 3: Top: contour plot of the uncorrelated proton spectral 
function of ^^O, Pmiss So(Pm.iss, Mmiss), sls Calculated with 
the Skyrme force and plotted vs. missing momentum Pmiss 
and missing mass, Mmiss = VV. Bottom: the correlated pro- 
ton spectral function of ^^O from [13], PmissSiiPmiss, Mmiss)- 
We show with black solid lines the cut Pmiss < 250 MeV and 
V < 640 MeV, quoted in Refs. [![ and 0] in both panels. 



shown in Figs. [3l |4] in contour plots; it has a similar 
behavior than the corresponding uncorrelated quantity 
but it is broader. One can see from Figs. [3l |4] that the 
correlated spectral functions in the two models exhibit 
differences; a detailed comparison of the two spectral 
functions is out of the scope of the present paper. Nev- 
ertheless, these differences do not affect our conclusions, 
as appears in the following, since most of the strength 
is concentrated along the stripe of maximum strength, 
even if the second model is more peaked at low momenta 
and it is more narrow around the center of the stripe. 
The center of the corresponding stripe obeys the follow- 
ing equation in the V, Pmiss plane for both models: 



V{P„ 



0.78 - 1.78 Pi, 
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where V is expressed in GeV^ and Pmiss in GeV/c, which 
is close to the approximate expression of Eq. ([I} . 

The concentration of the strength of the spectral func- 



FIG. 4: As in Fig. [31 but with the model spectral function of 
Ref. W^. 



tion in limited regions of space which project in some 
bands in the Mmiss, Pmiss plane suggests a complemen- 
tary analysis of the data specifically aimed at the decay 
of the ^^O protons. It consists in the following: to look 
for events which, in this plane, fall in one or several, de- 
pending of the accuracy of the data, regions of this plane 
selected to cover the lines of the maximum of the (un- 
correlated or correlated) strengths defined in Eqs. ^ [7]) 
so as to maximize the number of significant events while 
minimizing the background, i.e., the area. Correlated 
events can also be included in this way. The calculated 
proton spectral function in ^^O is sufficiently reliable, 
as it is established in connection with various nuclear 
physics experiments, to allow for this possibility. 

This kind of analysis, if feasible, precludes the subse- 
quent distortions of the pion kinematics after emission by 
the proton. We will comment later on that. A technical 
remark: for the correlated part, when transforming co- 
ordinates from {Pmiss , E) to {Pmiss , V) as in Eq. ([5]) , we 
impose that V stays positive, so we forbid a certain region 
of E, Pmiss space to be accessible; this in turn means the 
normalization of the correlated spectral function is not 
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FIG. 5: The total spectral function, So + ^i, within the two 
considered models. Top: the mixed model of So from the 
Skyrme force and Si from Ref. [l^; bottom: the model of 
Ref. fil. 



exactly 0.2 = / dEdPmiss Si{Prniss,E) but 



(8) 



We have checked that integrating over nega- 
tive values of V gives the missing normalization, 

J"^ dV J dPrmss Si{Prm,s,V) = 0.02. 

In Fig. [5] we show the total spectral function Sq + 
Si, in the two considered models, as a function of Pmiss 
and Mmiss, while in Fig. [5] we present the normalization 
integral of the spectral function 



dV / dPr, 



.SiPrmss,V). (9) 



It represents the number of events lost by applying a 
cut on the missing mass such that only the events which 
correspond to a missing mass larger than this particu- 
lar value Mmiss sue kept, irrespective of the momentum. 
If no cut on V is applied iV = 0.98. If we use the cut 
VV > 800 MeV, k < 250 MeV/c, presented in Figs. |3 
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FIG. 6: The normalization of the spectral function after the 
cut of Eq. m , as a function of the upper limit of integration 
in {Mmia3)max = ^/Vmax', the arrow shows the position of 
the proton mass, Mp. Black curves correspond to the model 
of Ref. and the red curves correspond to the model of 
Ref. [ll[; separate contributions are shown for the uncorre- 
lated/correlated part within both models. 



m and |5l the number of nuclear events is reduced by a 
factor of iV = 0.83 using the model of Ref. [Tol and by a 
factor of iV = 0.80 using the model of Ref. [ll| . However, 
future experiments are likely to have to introduce tighter 
cuts in order to reduce the background from the interac- 
tions of atmospheric neutrinos. If, for example, the cut 
VV > 900 MeV is imposed, ~ 44% of the events are re- 
moved, namely 26% (25%) from uncorrelated events and 
18% (19%), the near totality of correlated events in the 
first (second) considered model. With a tight constraint 
on Mmiss the fraction lost is quite appreciable if no other 
precaution is taken. The correlated analysis that we dis- 
cussed may allow a better efficiency. 

There are other effects which reduce the contribution 
of the bound nucleon decays. This includes a reduced 
phase volume which is (x ^/V for decays with production 
of light particles. Furthermore the very mechanism of the 
proton decay may be sensitive to the nuclear correlations. 
For example, if the decay amplitude is proportionalto the 
three quark wave function at the origin, see e.g. 



21|, the 



effect of suppression of the point-like configurations in 
bound nucleons [22, 23] would contribute, reducing the 
rate of the decay by about 14% for the Mp — VV = 100 
MeV cut. 

All the effects that we have discussed are genuine 
medium effects on the decay amplitude. They are not 
the whole story. The subsequent history of the pion, 
rescatterings / absorption in oxygen further reduces the 
number of "observable" pious. The inelastic scattering 
of pions clouds the message on the kinematics since the 
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inelastically scattered pion ejects a nucleon. The corre- 
sponding point in the M^iss j P-miss plane would be likely 
to fall outside the interesting regions delimited from the 
proton spectral function that we discussed in this work. 
Therefore for an analysis of the type suggested in this 
work, inelastically scattered pions may be considered as 
lost events. It may represent a reduction factor of about 
0.6. This is usually taken care of through a Monte Carlo 
evaluation which is beyond the scope of the present work. 

In conclusion we have introduced, in the problem of 
the identification of the decay events of the protons 
bound in the oxygen nucleus, the use of proton spectral 
function. It allows the prediction of the location in the 
Mmiss, Pmiss plane of the decay events. Our spectral 
function has an uncorrelated and a correlated part. It 
is has been tested against a number of nuclear physics 
experiments and the reliability of our prediction is 
sufhcient to be exploitable. We considered two models 
for the spectral function and the conclusions on the 



correlated cuts holds even if the two models exhibit 
some difference for the spectral hmctions. It appears 
that for the future nucleon decay experiments with tight 
cuts on the mass of the products of the proton decay 
it may be interesting to consider, as a complementary 
information for the decay events of the oxygen protons, 
correlated cuts on the mass and missing momentum 
obtained from the spectral function in order to decrease 
the background to signal ratio. The price to pay for this 
type of analysis is the loss of decay events where a pion 
produced in the decay is inelastically scattered, which 
clouds the reconstitution of the proton spectral function 
in oxygen. The loss in intensity however is moderate and 
it may be compensated by the advantage of a decrease 
of the background. 
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